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We examine a recently proposed nucleon-nucleon interaction, claimed by its authors both realistic 
and leading to a standing wave instability in symmetric nuclear matter. Contrary to these claims, 
we find that this interaction leads to a serious overbinding of *He, ^^O and ^''Ca nuclei when the 
Hartree-Fock method is properly applied. The resulting nuclear densities contradict the experimental 
data and all realistic Hartree-Fock results. 
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• Recently, a simple nucleon-nucleon interaction was proposed which is claimed both realistic and leading to a standing 

fSJ I wave instability in symmetric nuclear matter Although, strictly speaking, symmetric nuclear matter is a purely 
speculative object, it served for years as a testing ground for nuclear many-body theories and new insights into its 

' properties are of considerable interest. The hint that all these theories missed the spatial modulation of the nuclear 

^ matter density is provocative. When putting forward such claim, one has to make sure that the proposed interaction 

1 ' ■ satisfies constraints imposed by our knowledge of nuclear physics. 

[ The interaction considered in [ll| reads 
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y(ri, rs) = -aC{ri - r2)^e-'^'''-'''> + ^i^^<T>b[v^ - ra), (1) 
■ where < T > is the center-of-mass corrected average kinetic energy: 
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^ In the latter equation, single-particle orbitals (j>i relevant for the Hartree-Fock (HF) treatment are explicitly introduced. 
^ The auxiliary constant C = 27r^'^/^s^^/3, while the strength and range of attraction and the strength of contact 
' repulsion are chosen as a = 1690 McV fm^, s = 0.54 fm, f3 = 225 MeV 0. These parameters were intended to fit the 
' binding energy and the equilibrium density of nuclear matter and the binding energy of alpha particle (but see below). 
With these parameters, the authors reported reasonable values of the compressibility modulus of nuclear matter and 
of binding energies of even-even N = Z nuclei. 

It is crucial to understand that, although the authors refer to the HF method when describing their calculations for 
finite nuclei Q], they performed in fact only a very restricted minimization. This restriction is evident in a very small 
harmonic oscillator basis that has been used. In addition, the exchange integrals were not calculated, but assumed 
to be a fraction of the direct terms, depending on the average kinetic energy. Therefore, especially in view of quite 
important consequences claimed, an independent evaluation of the results of Ref. ||l| is called for. 

In this short note we report the results of regular HF calculations with the interaction (0) by which we determined 
binding energies of ^He, ^^O and '^''Ca nuclei. These binding energies are at variance with |l| and the obtained matter 
and charge densities of ^^O and '^''Ca are highly unusual. 



II. HARTREE-FOCK METHOD 

As the interaction (|l|) is spin- and isospin independent one assumes the fourfold degeneracy of the HF orbitals for 
even-even N = Z nuclei of interest. There are A/ A independent orbitals and we sum over them to obtain density 

p = '^flf I 4>i Neutron and proton densities are equal, Pn = Pp = 2p, and the total density equals 4p. Similarly, 
the total average kinetic energy < T > is quadruple of the sum of kinetic terms over independent orbitals. 
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The HF energy reads 

E^<T>+% j j d\id\2p{vi)p{Y2)Va{Yl - Tj) 
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-2^ / /rf'ridV20*(ri)</)*(r2)0j(ri),^,(r2)K(ri -r2) + 6/3 < T >^/^ j d^rp^ (3) 

where Va is the attractive part of (|l|). Remembering about differentiation of < T >, we obtain from d) a set of HF 
equations for the wave functions ipi and single-particle energies ei 

+ 3PV< T >p0, +aJ dV2p(r2)K(ri - r2)0^ 

j d\2<P*{v2)'Pj{vi)cp,{v2)Va{vl - V^) = Q0„ (4) 
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where the effective mass is given by m/m* = (1 — + 3(31 /{A\/ < T >)], with I = J d^rp^ . 

The numerical solution for (/i^ is straightforward, but tedious due to the exchange integrals. For spherically symmetric 
completely filled shells there is, however, a well known Slater method |^ to obtain the exact exchange potential. 
Adapting this general argument to the attractive potential Va we can express the exchange potential acting on the 
wave function 0„/ (r, 6*, 93) = Rn'i'{r)Yvm'{0,(p) as 

yaExi'n'l'm' — 

2/ + 1 f°° I 

E WTT^^^'^'^'^ / ^'^2rii?„i(^2)i?„'i'(r2)Vfe(ri,r2) (5) 
k=\i-i'\ + Jo J 

where the subshell index (nl) in the summation runs over the occupied orbitals. The coefficients A{k,l,l') are given 

by 

A{k, I, n) = y PkPiPn, (6) 

where Pi are Legendre polynomials. The functions Vk define the expansion of Va into spherical harmonics 

Fa(ri,r2) =^yfe(ri,r2)Ffe(cos0), (7) 
fe 

Vkirur^) = -aCs2e-('-?+'^i)/^'(2fc + l)[(fc + 1 + !i±^)/,(z) - zfk-i{z)l (8) 



where fk are the spherical Bessel functions of the imaginary argument, i.e. fk{z) — jk{iz) — \J j (2z)l}^j^x li{,z\ 
with /fc+i/2(2:) the modified Bessel function [|| and z ~ 2rir2/s^. 

III. RESULTS AND DISCUSSION 

In order to have a check on the HF results we used two different schemes, one spherical, using decomposition (|^), 
and the other three-dimensional. Both use wave functions defined over a spatial mesh. The general three-dimensional 
scheme being more time-consuming practically restricts the mesh size to about 30 points in each direction in one 
octant of space. The spherical scheme allows radial meshes of 100 points or more and may easily produce accurate 
solutions. 

We consider magic ^He, and ^"Ca nuclei for which the spherically symmetric solutions are expected. The HF 
problem was solved by the imaginary-time evolution. The convergence is rather slow for density and single particle 
energies, especially for ''°Ca. This is due to the buildup of central density peak which costs little energy in final stages 
of iteration. The three-dimensional scheme becomes impractical in this case but still its results tend to those of the 
spherical code. Below, we report densities calculted with the faster spherical code on the mesh of 100 points. 
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As the starting wave functions we took the results of [y. Therefore we could compare our initial energies and 
densities with those of [Q. We obtain for ^He the same energy as in [Q, but for ^^O and ^°Ca we find differences. 
These must be attributed to the error in energy introduced in ^ by the approximate treatment of the exchange 
integrals. Indeed, the integrals calculated analytically for s and p wave functions in -'^^O agree exactly with the results 
of our numerical codes. The correct values of the binding energy per nuclcon for initial configurations are 8.801 
MeV in ^^O and 11.137 MeV in ^^Ca, to be compared to 8.59 MeV and 10.76 MeV reported in 0. Thus, the exact 
calculation of the exchange integrals alone points to the overbinding problem with the interaction (|l|). This problem 
is magnified if one cares for the HF solutions. 

The first issue is the binding energy of alpha particle which bears on the determination of interaction constants 
The optimal wave function is more peaked than the gaussian used there. The HF binding energy per nucleon is 
8.076 MeV, i.e. 0.78 MeV more than in 0| where the experimental value corrected for Coulomb interaction was used. 
Thus, the interaction (|l|) overbinds already '*He when properly treated. 

The results for three nuclei are collected in Table 1. As seen there, the overbinding of ^®0, and especially '*°Ca, 
is really serious. For ^^O, the calculated binding is 173.57 MeV without Coulomb while the experimental value is 
127.619 MeV Allowing about 13 MeV for Coulomb energy (the direct term minus exchange, as it results from 
any realistic HF) we get more than 30 MeV overbinding. For ^"Ca, the calculated binding of 565.17 MeV, even after 
subtraction of about 71 MeV of Coulomb repulsion, is by about 152 MeV (!) larger than the experimental value of 
342.052 MeV In Table 1, we also give the difference in total binding (without Coulomb) between our results and 
those of 1^ to emphasise the importance of proper HF minimization. 

The calculated self-consistent nuclear densities are depicted in Figure 1. The tendency towards central peak 
development is evident. It is worth noting that this tendency is already seen in inaccurate results of [0. The density 
of the ^"Ca shown there exhibits a strange pile-up in the center. However, the HF results shown in Figure 1 allow to 
appreciate that this problem is even more grave: The central density is more than 1.5 times larger in the case of ^^O, 
and three times larger in ^"Ca than the experimental one (see e.g. pj). 



IV. CONCLUSIONS 

The exact HF calculations with the recently proposed interaction (^ for magic "^He, ^^O and ''°Ca nuclei show 
a serious overbindig problem. Associated nuclear densities develop central peaks taking a form unknown in nuclear 
physics. Both deficiencies grow with increasing mass. In view of the above results it is clear that the interaction 
proposed in |^] is very far from a realistic nucleon-nucleon force. Therefore, the assertions about the standing wave 
instability in nuclear matter made there, as related to unrealistic interaction, are unfounded. 
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TABLE 1 - Calculated HF binding energies per nucleon vs. results of ^ and the difference between total 

quantities (in [MeV]). 





B/A B/A in [] 


[\ B- B \] 


L| 


^He 

160 

40Ca 


8.08 7.3 3.11 
10.85 8.59 36.13 
14.13 10.76 134.77 



Figure caption 



Figure 1 

Total nuclear HF densities (thick lines) and densities from lit] (thin lines) for ^^O (dashed) and ^"Ca (solid). 
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